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Abstract. We study the Landau-Khalatnikov-Pradkin transformations 
(LKFT) in momentum space for the dynamically generated mass function 
in QED3. Starting from the Landau gauge results in the rainbow approxima- 
tion, we construct solutions in other covariant gauges. We confirm that the 
chiral condensate is gauge invariant as the structure of the LKFT predicts. 
We also check that the gauge dependence of the constituent fermion mass 
is considerably reduced as compared to the one obtained directly by solving 
SDE. 



1 Introduction 



Gauge theories of fundamental interactions have been highly successful in collat- 
ing experimental results in the perturbative regime. However, not all interesting 
phenomena can be understood in this approximation scheme. Confinement of 
quarks and gluons and the origin of hadronic masses are two examples. It is 
well known that on the distance scale of the order of hadron size, quarks behave 
as though their mass were 300 MeV, much larger than their vanishingly small 
current mass. In the context of covariant gauges, this effect may be attributed 
to an interplay of the behavior of gluon and ghost propagators in the infrared, 
[1]. These studies through Schwinger-Dyson equations (SDEs) have only been 
carried out in the Landau gauge. It is well known that if care is not taken in 
their truncation, solutions can be gauge dependent. 

Physically relevant strong dynamics is not restricted to QCD. There are also 
proposals of this kind to go beyond the standard model of particle physics such 
as (i) the technicolor models, (see[2] for an early review), (ii) top mode standard 
model and its modern variants ([3] and references therein) including those with 
extra dimensions, e.g. (5]. In view of these possibilities, it becomes all the 
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more important that the corresponding studies of SDEs be independent of the 
ambiguities stemming from the approximations employed. 

Let alone the complicated battle ground of non-abelian theories, the ambigu- 
ities associated with the lack of gauge invariance in the non-perturbative trun- 
cations of SDEs in abelian theories such as QED are also a challenging problem. 
Among other issues, it has ultraviolet divergences. Most of the numerical investi- 
gations of SDEs make use of the more convenient cut-off regularization method. 
However, such a cut-off can potentially introduce spurious gauge dependence [6], 
if employed naively. Consequently, one would hope that a neater testing ground 
for the validity of various truncation schemes could be provided by QED3, which 
is super-renormalizable. Nevertheless, a recent study, [7], has revealed that the 
problem persists even in this seemingly simple scenario. The main motivation of 
this work is to address this problem of gauge non- invariance in QED3. 

In the absence of fermion loops, i.e., in the quenched approximation, the lat- 
tice version of QED3 has been shown to be confining [8j, making its connection 
to the strongly interacting QCD. However, if the vacuum polarization tensor is 
taken into account at order a with massless fermions circulating in the loop, the 
potential is no longer confining, see Interestingly, for massive fermions, con- 
finement is reinstated, [TO] . These results have received numerical confirmation 

mm- 

In addition to being a toy model for SDE studies, QED3 is an attractive 
theory in its own right. In the field of superconductivity, see e.g. jl2j . it has been 
used in the context of high Tc superconductors. These studies have received a 
boost from new experiments [13]. QED3 has also found exciting applications in 
the recently explored unconventional quantum Hall effect in graphene [H]. Ad- 
ditionally, in the realm of dynamical generation of fundamental fermion masses, 
numerical findings on the lattice, results obtained by employing SDEs [15j and 
alternative methods [16], have yet to arrive at a final consensus, especially as 
regards the criticality of number of fiavors, and continue to provide a popular 
battle ground. 

The problem of gauge invariance in QED3 can be traced back to not em- 
ploying (or doing so incorrectly) the gauge identities such as the Ward-Green- 
Takahashi identities (WGTI) [17], the Nielsen identities [18] and the Landau- 
Khalatnikov-Fradkin transformations (LKFT) [T9]. In this article, we address 
this issue in the light of the LKFT. These transformations describe the specific 
manner in which Green functions transform under a variation of gauge. These 
are non-perturbative in nature, and are most readily defined in coordinate space. 
This was the basis for some earlier works [20] . Momentum space calculations are, 
however, more tedious [211 122[ [23 l 124] . Due to their non-perturbative nature, we 
expect LKFT not only to be satisfied at every order in perturbation theory, 
but also in phenomena which are realized non-perturbatively, such as dynamical 
chiral symmetry breaking. Initial steps were taken in [25] to apply these transfor- 
mations directly to the dynamically generated mass function. As a continuation 
of this effort, here we carry out an exact numerical exercise to study the behavior 
of the fermion propagator under these transformations in QED3. 

The traditional way to study gauge dependence is to make an ansatz for 
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the fermion-boson vertex, P ESI EZl EEl ES EQl EU [32] , see [33] for a brief re- 
view. The fermion propagator is then calculated by solving the resulting SDE 
in different covariant gauges. This method is useful to test how consistently a 
truncation scheme performs in different gauges. However, this line of attack does 
not naturally guarantee that LKFT for the fermion propagator will be satisfied 
gauge by gauge. In order to ensure that these transformations are satisfied, it 
is obligatory to study the constraints imposed by them on the dynamically gen- 
erated fermion propagator under a variation of gauge. One of the reasons it is 
important is because the said transformations guarantee that the condensate is 
a gauge-invariant quantity. Therefore they can provide a useful constraint on the 
truncation of SDEs. The way to do it is as follows. We start from the solution for 
the fermion propagator in the Landau gauge for a given truncation and simply 
perform a LKFT to find the result in any other gauge. The lesser the quanti- 
tative difference between the LKFT-generated fermion propagator and the one 
obtained by employing the truncation being tested in various gauges, the more 
reliable it is in terms of its gauge-covariant properties. 

We have organized this article as follows : In Sect. 2 we review the generalities 
behind the truncation of SDE in QED in general, and QED3 in particular. Sect. 3 
is devoted to the LKFT for the fermion propagator, and its behavior under 
gauge transformations as dictated by these transformations is reviewed in Sect. 4. 
Finally, we discuss our strategy and conclude in Sect. 5. 




Figure 1. SDEs for the fermion propagator, photon propagator and fermion-boson vertex 
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2 Truncation of Schwinger-Dyson equations 

We know that perturbation theory in QED does not generate fermion masses 
dynamically. Non-perturbative methods have to be employed and SDEs are the 
natural tool for calculation in continuum. The first three equations of the infinite 
tower of SDEs in QED are depicted in Fig. [H These correspond to the fermion 
propagator, the photon propagator and the fermion-boson vertex, respectively. 
The fermion propagator is coupled to the photon propagator and the fermion- 
boson vertex, which, through their own SDE, are coupled to the rest of the 
infinite tower containing higher point Green functions. 
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Wavefunction Renormalization in Various Gauges 
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Figure 2. Wavefunction Renormalization in the Rainbow Approximation in various gauges. 
The scale of the graph is set by = 1 



The SDE for the fermion propagator in QED3 is 



where q = k—p and is the dimensionful coupling of QED3. In this expression 
A^y{q) is the photon propagator and r^{k,p;^) the full fermion-boson vertex. 
The most general form of the (Euclidean) fermion propagator is 



F{p;0 



i^ + M{p;0 



(2) 
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where F is referred to as the fermion wavefunction renormalization and A4 is 
the mass function. The photon propagator in its general form can be written as 

. , . G{q^) ( , qtiqA , ^qi^qv _ ^ ^ , cip^qi^ (i\ 

A^u{q) = (^V ~ " '^"^^^ ^ ^ 

where Q is the photon wavefunction renormalization, which is gauge invariant, 
and ^ is the usual covariant gauge parameter. Expression ([1]) is a matrix equa- 
tion which can be converted into system of coupled non-linear scalar integral 
equations for M. and F after multiplying it, respectively, by 1 and ^ and taking 
trace. A favorite starting point in the quenched version of QED3, which con- 
sists in neglecting fermion loops (t/ = 1), is to choose a suitable ansatz for the 
fermion-boson vertex, in such a way that the SDE for the fermion propagator 
can be solved for the unknowns A4 and F. Possibly the simplest choice for the 
vertex, which allows us to understand the general features of dynamical chiral 
symmetry breaking, consists in replacing the fermion-boson vertex by its bare 
counterpart [9l[M]. This is the so-called rainbow approximation. 




To start with, the bare fermion propagator is considered massless. In this case 
the unknown functions defining the fermion propagator are found through the 
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following system of equations, 



1--^ / dk 



\ k'^ , 
1 ; m 

2kp 



F{p]i) 47r7o k-^ + M'^{k;i) 

M{p;0 _ «(e + 2) r ^^ kF{k-OM{k;0 



k + p 



■p 



k + p 



k — p 



(4) 



F(p;C) np Jo k^ + M^ik-0 

where a = e^/An as usual. In Figs. [2] and [3] we show the functions F and Ai 
as a function of momentum for various values of the gauge parameter. A4{p;£,) 
is roughly a constant for low momentum and falls as as p ^ oo [3l], for 

every value of ^, whereas F{p;(^) is a constant different from unity in the infrared 
in all but the Landau gauge, and approaches 1 at large p for all values of the 
gauge parameter. The non-trivial profile of the mass function is a result of chiral 
symmetry being dynamically broken. Observe that the curves corresponding to 
M{p',C) never cross each other for different values of ^. The resulting physical 
observables such as the chiral condensate (tpip) = — Tr5(x = 0;^), are gauge- 
dependent quantities, as can be seen in FigSl One can expect the source of such 
gauge dependence the fact that the bare vertex violates the WGTI, 

iqf,r^{k,p;0 = S-Hk;0 - S-Hp;0 ■ (5) 

In order to recover the reliability of a truncation to the SDE, one models the 
form of the three-point vertex, requiring, for example, the validity of the WGTI. 
Formally, in light of this identity, the vertex can be decomposed into two parts, 

r^{k,p-0=r;:{k,p;0 + r;^ik,p;0 , (6) 

such that qfj_rj^{k,p;^) = 0. rj^{k,p;^) remains undetermined by the WGTI 
and is called the transverse part. ^) is the longitudinal part. Although 

it is not uniquely fixed by the WGTI, a popular choice is the so-called Ball- 
Chiu vertex, Eq. (3.2) in [26]. An ansatz can be made for the transverse part of 
the vertex guided by perturbation theory, by the LKFT transformations, free- 
dom of kinematic singularities and its correct transformation properties under 
C, P and T. A few examples in this connection are the vertices proposed by 
Bashir, Curtis and Pennington, [27\ I29j. Yet the question remains, what should 
one expect to be the behavior of the fermion propagator in different gauges that 
ensures the gauge independence of physical observables associated with the phe- 
nomenon of DCSB? Whereas WGTI ensures the gauge invariance of the pole 
mass (in Minkowski space) , LKFT preserve the gauge independence of the chiral 
condensate. Therefore, in addition to the constraints of the WGTI, the ones of 
LKFT can also guide us toward improved truncations of SDEs. This is what we 
concretize in the next section. 

3 LKFT and the Fermion Propagator 

We start by writing the Euclidean space fermion propagator in both momentum 
and coordinate spaces in their most general forms : 
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Figure 4. Chiral condensate in the Rainbow approximation. The scale of the graph is given in 
units of lO'^e" = 1 

Expressions in Eq. d?]) are related through a d— dimensional Fourier transforma- 
tion. Assuming we know the fermion propagator in Landau gauge in momentum 
space, S{p;0), the LKFT relating the coordinate space fermion propagator in 
the Landau gauge to the one in an arbitrary covariant gauge reads 



This transformation leaves SDEs and WGTI functionally invariant in different 
gauges, and thus imposes strong constraints on the gauge covariance of the dy- 
namically generated fermion propagator. It is trivial to see that corresponding 
LKFT ensures the gauge invariance of the chiral condensate: 



S{x;C) 



S{x;0) ew[-i[Ad{0)-Ad{x)]] 



(8) 



where 




(9) 



TtS{x = 0; 

Tr5(x = 0; 0) exp [-i [AaiO) - A^ix = 0)]] 
TrS'(x = 0;0) = (V^V')o • 



(10) 
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Figure 5. Wavefunction Renormalization in the Rainbow Approximation in various gauges 
from LKFT. The scale of the graph is set by = 1 



Thus, the LKFT in SDEs studies should play an important role to address the 
issue of gauge invariance. Below we exemplify this considering QED3. 

The dynamically generated fermion propagator is usually known in momen- 
tum space. We can obtain a momentum-space representation for the LKFT in 
the following manner. Assuming we know S{p;0), we have the corresponding 
expressions for the propagator in coordinate space given by 



X{x;0) 
Y{x;0) 



d?k F{k;0) 
x2 J (27r)3 k^ + M^{k;0) 
(fk F{k;0) M{k;0) 



exp [—ik ■ x]k ■ x , 
exp [—ik ■ x] . 



11) 



(27r)3 k2 + M^ik;0) 
Performing the angular integrations, these expressions get simplified as follows. 



X{x;0) 
Y{x;0) 



1 



27r2x2 70 
1 



k2+M^{k;0) 
2F(fc;0) M{k;0) 



cos kx 



sin kx 



kx 



27r2 7o k^ + M^{k-{)) 



sin kx 

kx 



(12) 



It is straightforward to check that for d = 3, /^3(0) — A^{x) = —ia^/2. Hence, 
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S{x;0 = 5(x;0)e-"^ (13) 

with a = After applying LKFT to the expressions in Eq. (I12p and Fourier 

transforming the results back to momentum space, we get 

- - f d^x exp [ip ■ x]X{x- e) , 
P J 

-- - J cfx exp[ip ■ x]Y{x;^) . (14) 



p2 _|__yV42(^p.^) 
F{p;OM{p;0 

pi + M'^{p]i) 



Mass Function from LKFT in Various Gauges 
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Figure 6. Mass Function in the Rainbow Approximation in various gauges from LKFT. The 
scale of the graph is set by = 1 



Angular integrations lead to the final momentum space representation of the 
LKFT for the fermion propagator (provided a > 0) : 



F{p-i) 



p2 _^ _/V42(p;^) 

F{p-0 M{p;0 
p2 M'^{p;i) 



dk k 



F{k;0) 



a 

T^P Jo 



— dkk 



k2+M^ik;0) 

F{k;0) M{k;0) 
k^+M^k;0) 



1 1 1 



A- 



A+ 



1 



1 

A+ 



(15) 
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where = + (/cztp)^. Thus knowledge of S{p; 0) is the input required to ob- 
tain the same in an arbitrary covariant gauge. One can further confirm that this 
representation preserves the gauge invariance of the chiral condensate upon inte- 
gration the second expression over J d?p/ . The left-hand side corresponds 
to the integrand of the condensate in momentum space. 

4 Behavior of the fermion propagator 

We shall now review the behavior of the fermion propagator in different gauges 
as dictated by the corresponding LKFT. Starting with the chirally asymmetric 
solution of the fermion propagator in the Landau gauge, we apply the LKFT to 
obtain the off-shell fermion propagator in an arbitrary covariant gauge. 

Condensate Integrand in Various Gauges 

0.05 I ' ^ ' -^T ^ -^T ^ -^T ^ ^ ' ^ 




0.001 0.01 0.1 1 10 100 1000 

p 

Figure 7. Condensate integrand in the Rainbow Approximation in various gauges from LKFT. 
The scale of the graph is set by = 1 

The behavior of the fermion propagator is depicted in Figs [5] and [6l Observe that 
the wavefunction renormalization in qualitatively the same both for the SDE 
solution in different gauges and the LKF-transformed solution. For the case of 
the mass function, the asymptotic behavior remains qualitatively the same, as the 
analytical treatment in Ref. [25] confirms, but the key observation is something 
that might have been anticipated but has not before been demonstrated; namely, 
a rearrangement in the mass function in different gauges in such a fashion as to 
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render the chiral condensate gauge invariant. This is illustrated in Fig. [TJ where 
we depicted the condensate integrand 

as a function of momentum for various values of the gauge parameter. A simple 
numerical exercise reveals that the area under each curve is equal, rendering the 
condensate gauge invariant. This can explicitly be seen in Fig. [8l where the flat 
line corresponds to the gauge (in)dependence of the condensate calculated from 
the LKF transformed solutions as compared with the SDE results. 




Figure 8. Gauge dependence of the condensate in the Rainbow Approximation from SDE and 
LKFT. The scale of the graph is set by = 1 



Constraints on the behavior of the fermion propagator arising from LKFT 
have implications on other physical observables, like the constituent mass in 
Euclidean space defined hy m = A4{p = m). A comparison from the results ob- 
tained from the rainbow truncation of SDE [31] and those coming from LKFT is 
depicted in Fig. 9. Gauge dependence of the constituent mass arising from LKFT 
has been reduced considerably, as compared with the sharp gauge dependence 
of the SDE result. 

It has been shown that the fermion propagator has complex mass like singu- 
larities implying confinement ^Jj . Through the evaluation of the Euclidean-time 
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Constituent Mass in Various Gauges 
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Figure 9. Constituent mass in the Rainbow Approximation in various gauges from LKFT. The 
scale of the graph is set by = 1 



Schwinger function, one can determine whether or not the fermion propagator 
corresponds to a physically observable state or not [Ml [36] . It is defined as 



(fx 



exp [i(pot + P • x)] a{p ;0 , 



where 



(27r)3 



(17) 



(18) 



p"^ + M'^{p]i) 

It can be shown that if there is a stable asymptotic state of mass m associated 
with this propagator, then 

for large Euclidean t. This readily implies 



— lim ln[Z\(t,^)] = m . 
t^oo at 



(19) 



On the other hand, two complex conjugate mass like singularities with complex 
masses n = a±ib lead to an oscillating behavior such as 

Z\(t;0 ^e~"* cos(5t + 5) (20) 
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for large t. If we want to associate physical interpretation to these masses, there 
gauge independence needs to be established. As the above-mentioned method is 
not very accurate, associating exact numbers to a and b in various gauges might 
not reflect the gauge (in) dependence of these quantities. However, we present a 
qualitative idea through Fig. 10. 



eg 



Confinement Test in Various Gauges 



-10 



LKF vs. SDE 



T 1 1 

SDE ^=0 - 
SDE 1=1 o 
LKFT 1=1 



10 



?A A A AAA A A 
s aA o 

O 

O 

O 

O 



15 



20 



25 



30 



35 



40 



45 



50 



t 



Figure 10. ln[A{t; ^)/2] for C = 0, 1 through SDE and LKFT generated mass function 



We see that the LKFT-generated ln[Z\(t; l)/2] is qualitatively much closer to 
ln[Z\(t; 0)] than the SDE-generated ln[Z\(t; l)]/2. Exact numbers require solving 
SDE for complex momenta. 

5 Discussion and Conclusions 

To conclude, just as WGTI provides useful constraints on the truncations of 
SDEs, so do the LKF transformations. It is constructive to compare the solutions 
of the SDE in different gauges against those generated by the LKFT, p5] . If the 
off-shell propagators obtained through these two methods differ drastically from 
each other, the truncation adopted does not conform to the gauge covariance 
properties of the fermion propagator. For example, in Fig. [TT] we compare the 
gauge dependence of the condensate for two different truncation schemes with 
our strategy. On the one hand we show the rainbow approximation, where results 
from SDE and LKFT show a clear difference. On the other hand, we show the 
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same comparison for the Ball-Chiu vertex [26]. This proposal fulfihs the WGTI, 
hence rendering a less pronounced gauge dependence for the condensate, as can 
be seen from the small difference between SDE and LKFT results. This procedure 
detailed in the paper is also valid when one goes beyond the quenched truncation 
and tries to calculate the vacuum polarization. The same vertex must appear 
in both the fermion and photon gap equations. The bare vertex is inadequate 
because the vacuum polarization tensor must be transverse. However, once a 
vertex is constructed to ensure this, LKFT will provide the vertex required in 
another gauge. Therefore, more definitive conclusions upon the most appropriate 
vertex ansatz can be drawn after a detailed study of the LKFT for the vertex 
function itself [37]. This work is in progress. A natural extension of this work is 
to study the gauge covariance of the quark propagator through the generalized 
LKFT of QCD. 
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Figure 11. Comparison of the gauge dependence of the condensate for the Rainbow Approxi- 
mation and the Ball-Chiu vertex from SDE and LKFT. The scale of the graph is set by = 1 
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